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ABSTRACT 


The  irregular  space-time  sampling  of  any  finite  region  by  an  orbiting  satellite  raises  difficult  questions  as  to 
which  frequencies  and  wavenumbers  can  be  determined  and  which  will  alias  into  others.  Conventional  sampling 
theorems  must  be  extended  to  account  for  both  irregular  data  distributions  and  observational  noise — the  sampling 
irregularity  making  the  system  much  more  susceptible  to  noise  than  in  regularly  sampled  cases.  The  problem 
is  formulated  hoe  in  terms  erf' least-squares  and  applied  to  spacecraft  in  10-day  and  17-day  repeating  orbits. 
The  "mamond^pattem  "laid  down  spatially  in  such  repeating  orbits  means  that  either  repeat  period  adequately 
samples  the  spatial  variables,  but  the  slow  overall  temporal  coverage  in  the  1 7-day  pattern  leads  to  much  greater 
uncertainty  than  in  the  shorter  repeat  cycle.  The  result  is  not  definitive  and  it  is  not  concluded  that  a  10-day 
orbit  repeat  is  the  most  appropriate  one.  A  major  conclusion  however,  is  that  different  orbital  choices  have 
potentially  quite  different  sampling  characteristics  which  need  to  be  analyzed  in  terms  of  the  spectral  characteristics 
of  the  moving  sea  surface.  Conclusions  drawn  from  model  asrimilatiop  studies  need  to  be  placed  in  a  context 
of  the  reality  of  their  spectral  content  via4-vis  the  ocean.  r  TS  i  C  thdL  I  ^ — 

A/aao/v-esr-  T-/s y / 


of  all  cases.  Consider  an  even  number,  N,  of  obser¬ 
vations  f(nAq),  n  =  0  to  N  -  1  spaced  equally  in  a 
continuous  coordinate  q,  at  intervals  Aq  (q  is  delib¬ 
erately  made  abstract;  it  can  be  space  or  time).  Then 
(e.g.,  Hamming  1973),  the  expansion  functions 


1. 

Subsatellite  tracks  lay  out  a  geographically  and  tem¬ 
porally  complex  sampling  pattern  on  the  sea  surface. 
As  with  any  measurement  system,  understanding  of 
its  capabilities  and  discussion  of  orbital  choice  neces¬ 
sitates  at  least  a  semiquantitati  ve  determination  of  its 
sampling  characteristics.  For  an  instrument  making 
point  measurements,  uniformly  spaced  in  time,  one 
normally  determines  the  Nyquist  frequency  from  the 
known  characteristics  of  the  physical  process  being 
measured,  and  the  physical  characteristics  of  the  mea¬ 
surement  system.  Usually  there  is  some  alias-band  of 
frequencies  which  one  cannot  hope  to  observe  properly, 
and  one  designs  the  system  with  that  in  mind.  A  similar 
determination — finding  which  frequencies  and  wave- 
numbers  are  determinable,  and  which  are  aliased,  from 
a  given  orbit — is  an  essential  first  step  before  one  can 
discuss  alternative  choices.  The  immediate  motivation 
is  the  need  specifically  to  understand  altimetric  sam¬ 
pling,  but  the  results  apply  equally  well  to  scatterom- 
eter,  or  infrared,  or  in  general,  to  any  other  instru¬ 
ment  sampling  three  dimensions  in  a  nonregular 
manner. 

Most  conventional  sampling  theorems  (e.g..  Free¬ 
man  1963 )  are  derived  for  uniformly  spaced  data.  Al¬ 
though  the  satellite  case  is  much  more  complex,  it  is 
helpful  to  recapitulate  the  most  basic  and  important 
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where  L  =  NAq,  are  a  complete  basis  for  the  samples: 


ao 


"n  2xr(nAq) 


f(nAq)  =  —  +  2  a.- cos 


.  <w4>-'  L  .  2 xr(nAq) 
+  2  b,  sin - - - 

T- 1  ^ 


n  =  0  to  N  -  1  (2a) 


exactly,  if 


N- 1 


Nak  =  2  /(«A$)cos|y  knAq'j  , 

N~ 1  /  2x  \ 

Nbk  =  2  f{nAq)  sin^  —  knAq j  .  (2b) 
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(The  normalization  is  2 N for  k  =  0,  N. )  Any  frequencies, 
u),  outside  the  “baseband,”  |w|  2*/2Aq,  alias  into 

the  set  of  frequencies. 
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involved  id  (1). 

>  This  “ShaAimn- Whittaker  sampling  theorem”  is 
sometimes  best  ,vie wed  from  a  converse  point  of  view: 
{h&t  if  A  0  for  |<o|  >  (2x/2A q),  (a  “band- 

limited”  functiojk.)  then  if  there  are  an  infinite  number 
otabservatidQv/C?)  can  be  reconstructed  perfectly 
from  the  interpolation  formula 


f(Q)=  2  f(nAq) 


(4) 


If  the  band-limiting  is  violated,  one  can  still  use  the 
formula  (4),  but  the  reconstructed  continuous  function 
f(q)  does  not  coincide  with  the  correct  function  except 
at  the  observation  points  themselves.  High  frequencies 
alias  (i.e.,  appear  as  lower  frequencies)  with  potentially 
disastrous  consequences  when  one  goes  to  interpret  the 
data,  as  depicted  in  Fig.  1 . 

The  derivations  of  the  aliasing  formula  and  the  in¬ 
terpolation  rule  (4)  are  elementary  for  equally  spaced 
data  (e.g.,  Bracewell  1978).  But  subsatellite  data  lo¬ 
cations  while  not  randomly  distributed,  are  quite  ir¬ 
regular,  and  the  normal  derivations  will  not  apply. 
Furthermore,  the  presence  of  noise  in  the  data  is  not 
easily  accommodated  by  the  conventional  result.  We 


need  a  different  approach,  which  we  will  try  to  elucidate 
by  first  applying  it  to  the  known  case  just  described. 


2.  Sampling  theorems  by  least  squares 

If  we  attempt  a  least-squares  fit  of  the  basis  functions  - 
( 1 )  to  f(q)  in  the  form 

(2irk 


mm: 


f(nAq)  -  2  ak  cos  —  nAq 


\  L 


-  2  bu  sin 

k 


(5) 


it  is  well  known  that  the  conventional  Fourier  coeffi¬ 
cients  are  obtained.  The  least-squares  approach  would 
not  normally  be  used  for  the  equally  spaced  case,  be¬ 
cause  it  fails  to  take  advantage  of  the  enormous  com¬ 
putational  efficiency  made  possible  by  the  orthoi- 
gonality  of  the  basis  functions  and  fast  transform 
techniques.  But  least-squares  has  the  advantage  of  gen¬ 
eralizing  to  the  cases  we  need  to  examine. 

Even  though  the  fit  (5)  is  exact,  with  a  minimum 
of  zero  if  the  full  basis  set  ( 1 )  is  used,  and  ordinary 
numerical  Fourier  transforms  do  a  fine  job  on  finding 
the  ak,  bk,  let  us  rewrite  (5)  as:  Find  the  best  (in  the 
minimum  mean-square  residual  sense)  solution  to  Eq. 
(2a),  which  we  rewrite  as 


Aa  =  0, 


(6) 


where 
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Qni  2 
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,  0  = 

f[2Aq] 

f[(N-  1  )Aq\ 

We  solve  (6)  using  the  singular  value  decomposition 
(SVD);  for  example,  see  Lawson  and  Hanson  ( 1974) 
or  Wunsch  ( 1978)  and  write  the  solution  as 

*  uT8  N 

a  =  2  x  v<+  ^  (7) 

<-i  <-*+i 

with  N,  where  u,  and  y,  are  the  column  vectors  of 


the  singular  value  decomposition  of  A  making  up  the 
columns  of  U  and  V  respectively,  where 

A  =  UAV  r.  (8) 

The  X,  are  the  “singular  values,”  the  elements  of  the 
diagonal  matrix  A .  We  write  a  to  distinguish  the  com¬ 
puted  solution  from  the  true  one.  The  y,  are  constants. 

[  If  the  summations  over  a*  and  bk  are  carried  to  the 
upper  limits  of  ( 2a )  the  problem  as  posed  in  ( 6 )  is  fully 
determined  and  the  SVD  machinery  is  unnecessary. 
One  of  its  great  advantages  however,  is  that  this  con¬ 
ventional  situation  is  included  as  a  special  case  of  the 
much  more  general  problems  encountered  below.]  We 
adopt  the  usual  convention  that  rows  and  columns  of 
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FIG.  1.  The  classical  aliasing  effect  of  uniform  sampling  too  slow 
to  distinguish  a  frequency  above  the  Nyquist  frequency  from  one 
below,  as  in  expression  ( 1 ).  The  interpolated  function  agrees  with 
the  original  function  at  the  sample  points,  and  is  smooth  in  between. 
But  any  estimate  of  the  derivatives  of  the  original  function  made 
from  such  a  grossly  undersampled  situation  will  be  erroneous. 

the  diagonal  matrix  A  which  are  all  zeroes  are  dropped, 
as  are  the  corresponding  columns  of  the  U  and  V  ma¬ 
trices  [i.e.,  the  null-space,  if  present,  is  dropped  from 
the  representation  (8)].  Expression  (7)  is  completely 
general,  the  second  terms  on  the  right  representing  the 
null-space:  the  components  with  unknown  and  inde¬ 
terminate  coefficients  7,.  The  present  case  is  special 
because  the  singular  values  X,  are  all  about  the  same 
size,  and  there  are  k  =  N  of  them,  just  the  number 
necessary  to  solve  the  problem.  Figure  2a  displays  the 
singular  values  of  the  present  case. 

A  part  of  the  SVD  machinery  is  the  so-called  reso¬ 
lution  matrix  W  r.  We  cannot  provide  a  full  discussion 
of  the  interpretation  of  this  operator  here  (see  Wiggins 
1972;  Wunsch  1978 ),  but  the  following  may  help  those 
unfamiliar  with  it. 

If  the  null-space  of  the  solution  vanishes,  then  all 
elements  of  the  solution  are  completely  determined 
and  denoted  “fully  resolved.”  If  a  null-space  exists,  it 
is  a  straightforward  matter  to  show  that  the  relationship 
between  the  correct  solution  a  and  the  one  determined, 
«,  is 

a  =  W  Ta.  (9) 

Equation  (9)  is  readily  interpreted  as  showing  that  in 
the  general  case,  the  solution  estimated  is  a  linear  com¬ 
bination  (an  average)  of  the  elements  of  the  true  so¬ 
lution.  If  element  /  of  the  true  solution  is  fully  resolved, 
then  the  corresponding  row  i  of  W  T  is  the  ith  row 
of  the  identity  matrix.  It  is  an  elementary  demonstra¬ 
tion  that  in  the  fully  resolved  case,  the  entire  W 7  ma¬ 
trix  reduces  to  the  identity,  as  it  must.  In  the  opposite 
extreme,  when  the  data  contain  no  information  about 
element  1  of  the  solution,  the  z'th  row  of  VV  T  is  all 
zeroes. 


In  the  present  case  everything  is  fully  determined, 
and  W  T  is  the  identity.  The  problem  was  solved 
slightly  differently  from  the  way  we  have  posed  it, 
though.  In  particular,  because  the  computer  coding  was 
slightly  simpler,  and  because  it  is  illustrative  of  a  more 
general  machinery,  the  basis  functions  employed  in¬ 
cluded  the  sine  of  both  zero  frequency,  and  the  Nyquist 
frequency — neither  of  which  belongs  to  the  basis  set 
( 1 ).  Figure  2b  displays  the  diagonal  elements  of  W  r 


Fig.  2.  (a)  Singular  values  of  the  fit  of  16  points  to  their  Fourier 
series  by  least-squares.  Two  of  the  values  are  slightly  higher  than  the 
others  because  the  mean  square  value  of  unity  and  the  cosine  of  the 
Nyquist  frequency  are  2X  those  of  a  cosine  or  sine.  This  phenomenon 
could  be  removed  by  appropriate  column  scaling,  but  it  was  not 
thought  worthwhile  to  do  so.  (b)  Diagonal  component  of  the  reso¬ 
lution  matrix  for  the  same  fit  as  in  Panel  (a).  All  values  are  unity 
( full  resolution )  except  for  two  that  vanish.  These  latter  two  are  an 
artifact  deriving  from  carrying  two  extra,  useless  functions  in  the  fit. 
notably  the  sine  of  zero  frequency  and  the  sine  of  the  Nyquist  fre¬ 
quency  [neither  belongs  to  the  set  in  Eq.  (I)].  One  of  the  great 
virtues  of  the  SVD  is  its  ability  to  automatically  handle  this  indeter¬ 
minancy  (we  had  16  equations  in  18  unknowns)  which  was  done 
solely  because  the  computer  coding  was  marginally  simpler.  The 
handling  is  automatic  in  that  one  is  told  that  no  information  is  avail¬ 
able  about  these  two  coefficients,  and  if  one  solved  for  them,  the 
SVD  would  produce  a  best  estimate  of  zero. 
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for  the  case  actually  run;  one  sees  that  all  of  them  are 
unity  (the  remaining  elements  in  each  of  the  corre¬ 
sponding  rows  are  zeroes),  with  the  exception  of  the 
two  corresponding  to  the  two  extra  sinusoids — for 
which  the  elements  are  zeroes — showing  that  the  data 
distribution  we  had  contains  no  information  whatever 
about  these  two  basis  elements.  An  important  property 
of  the  SVD  method  is  its  ability  to  carry  two  unnec¬ 
essary  elements  and  to  inform  the  user  that  they  are 
unnecessary.  Later  on,  we  will  see  more  interesting 
cases. 

All  least-squares  solutions  contain  a  priori  assump¬ 
tions  about  noise  and  solution  statistics  ( Wunsch  1989) 
and  the  particular  ones  leading  to  (7)  are:  uniform 


b  FREQUENCY  RADIANS/UNIT  TIME 

FIG.  3.  (a)  Diagonal  of  the  resolution  matrix  of  regularly  spaced, 
aliased  case.  Apart  from  the  cosine  of  frequency  *7  A q  all  nonzero 
resolutions  have  been  reduced  to  0.5.  (b)  One  of  the  rows  of  the 
resolution  matrix  corresponding  to  3b,  showing  that  the  cosine  coef¬ 
ficient  at  that  particular  frequency  is  determined  only  in  linear  com¬ 
bination  with  the  cosine  of  the  alias  frequency  { elements  correspond¬ 
ing  to  the  sine  (arts  have  been  suppressed;  they  all  vanish).  The 
result  is  an  inability  to  distinguish  the  frequency  from  its  alias  ( without 
provision  of  further  information  I . 


variance  in  the  solution  elements,  and  zero  covariance 
between  solution  elements;  also,  any  noise  in  the  ob¬ 
servation  is  uncorrelated  from  one  measurement  to 
another,  and  is  of  constant  variance. 

It  is  instructive,  before  moving  on,  to  examine  mod¬ 
ifications  of  the  present  example. 

(i)  Suppose  we  try  to  resolve  frequencies  outside 
the  fundamental  baseband  |a>|  (2*/2A q)l  That  is, 
suppose  we  seek  to  find  the  Fourier  coefficients  out  to 
frequencies  =  2  X  Nyquist  frequency  =  ( 2ir/L )  N! 

The  problem  is  now  underdetermined,  because  we 
are  looking  for  twice  as  many  Fourier  coefficients  as 
data  points,  but  as  will  be  seen,  underdetermination 
has  to  be  faced  later  anyway.  The  diagonal  of  the  res¬ 
olution  matrix  W  T  for  this  case  is  displayed  in  Fig. 
3a  and  a  typical  row  in  (3b).  The  resolution  of  the 
sines  of  frequencies  0,-ir/Aq  and  2*7  A<?  is  zero  as  ex¬ 
pected.  The  resolution  of  the  amplitude  of  the  cosine 
of  frequency  *7  A q  remains  unity,  but  none  of  the  other 
coefficients  is  now  fully  resolved.  The  diagonal  element, 
which  was  unity  before,  is  now  reduced  to  'A,  and  'A 
appears  in  the  row  at  corresponding  frequency 
+  (2*7 Atf)  which  is  precisely  the  alias  frequency  (3). 
The  form  (9)  thus  states  that  the  information  available 
is  only  adequate  to  determine  the  average  of  the  ;th 
Fourier  coefficient  in  the  baseband  and  that  of  its  higher 
frequency  alias  in  the  form  o,  =  |(«,  +  «aiias0r,)-  If  the 
baseband  frequency  were  truly  absent  and  the  higher 
frequency  present,  then  analysis  would  nonetheless 
produce  a  nonzero  value  within  the  baseband.  In  par¬ 
ticular,  the  coefficient  of  the  cosine  of  zero  frequency, 
the  mean,  is  determined  as  the  sum  of  the  true  coef¬ 
ficient  of  that  frequency  and  that  at  frequency  2*/  A q — 
one  of  its  aliases.  The  rows  of  W  T  thus  produce  an 
explicit  statement  of  all  the  frequencies  contributing 
to  the  estimates  of  Fourier  amplitudes  at  any  given 
frequency.  In  the  present  uniformly  distributed  noise- 
free  case,  the  aliases  are  simple  to  understand;  when 
we  get  to  irregular  noisy  data,  the  dependences  become 
very  complex. 

Practitioners  of  inverse  methods  will  recognize  that 
we  could  give  a  larger  a  priori  weight  to  the  lower  fre¬ 
quencies  if  we  chose,  thus  assigning  all  of  the  amplitude 
to  it.  Assignment  of  such  a  priori  weights  is  outside  the 
normal  scope  of  sampling  theorems,  however,  which 
treat  all  frequencies  on  an  equal  footing. 

(ii)  Now  revert  to  consideration  of  the  baseband 
alone,  but  take  one  of  the  original  uniformly  spaced 
samples,  and  remove  it,  replacing  it  by  a  sample  at  an 
irregular  interval  (see  Fig.  4).  Freeman  ( 1965)  shows 
how  to  prove  that  the  irregularly  spaced  sample  is  a 
complete  substitute  for  the  missing  regular  one.  But  he 
deals  only  with  the  case  of  perfect  observations,  and  1 
wish  to  demonstrate  the  effect  of  noise. 

We  replace  one  sample  f(mAq)  by  a  sample  at  an 
irregularly  spaced  point /(/m),  and  the  corresponding 
missing  equation  in  (2a)  or  (6)  is  replaced  by 
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Flo .  4.  Point  number  7  has  been  displaced  from  its  “regu¬ 
larly  spaced”  location  to  a  position  much  closer  to  its  neighboring 

pt.-.:»  6. 


„  ,  2*r 

fOm)  =  -T  +  2  arCOS  —  t„ 

1  r=l  L 


<A'/2)-l 
+  2 


r«  1 


,  .  2irr 

b,  sin  — 1„ 


Formally,  there  are  still  TV  equations  in  TV  unknowns 
and  this  observation  is  the  basis  of  the  conclusion  (for¬ 
malized  by  Freeman  1965),  that  one  does  not  need 
regularly  spaced  data  to  preserve  the  perfect  recon¬ 
struction  in  (4).  We  persist  however,  in  solving  it  by 
the  S  VD;  the  solution  is  still  then  of  the  form  ( 7 ) .  Use 
of  the  SVD  demonstrates  the  real  effect  of  the  irregular 
spacing:  one  of  the  singular  values  (displayed  in  Fig. 
5)  has  become  very  small. 

Because  the  covariance  of  the  solution  is 

<(«  —  <«>)(«  ~  <“»r)  =  <r2VA~2Vr,  (10) 

where  a2  is  the  observation  noise  variance  ( for  example, 
see  Wunsch  1978),  the  effect  of  the  irregularly  spaced 
observation  is  to  reduce  the  effective  rank  in  the  pres¬ 
ence  of  noise.  As  the  spacing  becomes  more  irregular, 
this  noise  susceptibility  grows.  If  the  effective  rank  is 
significantly  reduced,  some  of  the  Fourier  coefficients 
will  become  indeterminate,  or  effective  aliases  of  others, 
or  simultaneously  both  of  these  things. 

As  the  number  of  irrregularly  spaced  points  in¬ 
creases,  and  as  the  spacing  between  any  pair  of  points 
is  reduced,  the  singular  values  decline.  Thus  the  limiting 
case,  when  all  the  samples  are  squeezed  arbitrarily  close 
together  in  an  arbitrarily  small  interval  remains  a  hy¬ 
pothetical  sampling  strategy  for  perfect  data,  but  is 
ruled  out  in  practice  because  the  slightest  observational 
noise  would  render  arbitrarily  large  the  variance  of  the 
result  and  the  numerical  values  as  useless.  (The  limit 
is  tricky  in  any  case;  see  Jerri  1977.) 

When  we  move  to  two  and  three  dimensions,  the 
SVD  remains  the  analysis  tool  of  choice  because  in  a 


precise  sense  it  provides  all  possible  information  con¬ 
cerning  the  solution  (for  example,  a  review  of  its  prop¬ 
erties  may  be  found  in  Wunsch  1978).  But  the  com¬ 
putational  load  grows  very  rapidly,  and  the  amount  of 
information  is  so  large  as  to  be  overwhelming,  and 
indeed  more  than  we  need  for  present  purposes. 

3.  Tapered  least-squares  approach 

In  any  estimation  problem,  there  are  generally  two 
types  of  errors:  (i)  those  owing  to  wholly  missing  in¬ 
formation,  and  (ii),  those  owing  to  noise  in  the  obser¬ 
vations  (a  third  error,  often  the  most  important,  is 
model  failure;  but  it  is  not  an  issue  here ) .  It  is  ( i )  which 
appeared  in  section  2,  when  we  sought  twice  as  many 
Fourier  components  as  we  were  entitled  to  estimate 
with  the  available  data,  and  is  the  error  owing  to  a 
“failure  to  resolve.”  When  computed  by  an  SVD,  the 
unresolved  components  are  assigned  zero  amplitude; 
i.e.,  the  null  space  coefficients  are  set  to  zero.  If  the 
measurements  are  not  perfect,  the  noise  gives  rise  to 
an  uncertainty  in  those  components  which  have  been 
determined.  This  latter  noise  is  what  appears  as  the 
solution  uncertainty  in  conventional  overdetermined 
least-squares.  The  advantage  of  the  SVD  is  that  it  is 
able  to  compute  these  errors  separately  (and  because 
their  origins  and  cures  are  quite  different,  a  separation 
is  often  extremely  helpful);  the  main  disadvantage  is 
the  large  computational  load  involved. 

Conventional  “tapered  least-squares”  combines 
these  errors  into  a  single  error  estimate  which  is  more 
easily  affordable.  The  result  is  obtained  most  simply 
from  the  Gauss-Markov  theorem,  in  the  restricted 
form  in  sections  5-7  of  Liebelt  ( 1967).  To  use  this 
result,  we  must  however,  specify  a  prior  estimate  of 


Fig.  5.  Singular  values  of  Fourier  fit  to  data  distributed  as  in  Fig. 
4;  singular  values  change  from  those  displayed  in  Fig.  2a.  with  one 
of  them  becoming  very  small — thus  making  the  result  vulnerable  to 
noise. 
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the  variance  of  the  solution  to  (6)  relative  to  the  noise 
in  the  observations. 

Suppose  that  the  observational  noise  in  the  ith  ob¬ 
servation  is  rtj  and  has  an  “a  priori  covariance”: 

{riiTij)  =  <r%, 

with  (a,aj)  =  A2S0  as  the  a  priori  covariances  of  the 
Fourier  coefficients;  then  we  have 


a  =  |atA  +  I  j  V/? 

(Ha) 

r  a2  I-1 

e-„’[a'a  +  pi]  , 

(lib) 

o=Ar  AAr+^I  0 

(12a) 

V  —  1 

E  =  A2  I  —  Ar|  AAr  +  I 

)  A  ,  (12b) 

where  E  =  ((a  -  ( a))(a  -  (a))T). 

Expressions  (11)  and  ( 12)  are  algebraically  identical; 
a  choice  between  them  is  made  for  numerical  efficiency 
depending  upon  whether  there  are  more  observations 
than  unknowns  [use  (11)],  or  more  unknowns  than 
observations  [use  ( 12)] — choosing  so  as  to  minimize 
the  dimension  of  the  matrix  to  be  inverted. 

To  demonstrate  how  E  combines  the  overall  mean- 
square  errors,  we  display  in  Fig.  6a  the  result  for  the 
one  irregularly  spaced  point  in  Fig.  2,  but  where  now, 
we  have  asserted  that  (a2/ A2)  =  0.1.  Figure  6  shows 
the  normalized  variance  of  the  cosine  coefficients  for 
the  classical  case  (the  covariance  matrix  is  diagonal). 
For  the  case  of  one  (now  noisy)  irregularly  spaced  ob¬ 
servation,  the  variance  increases  as  shown  in  the  same 
figure.  There  are  off-diagonal  components  in  this  co- 
variance  matrix  (Fig.  6b).  For  this  particular  pair  of 
frequencies,  there  is  a  negative  correlation  between  the 
estimation  errors  of  each.  As  the  number  of  data  spac¬ 
ing  irregularities  increases,  the  variance  generally  grows 
monotonicaily,  with  some  wavenumbers  being  more 
sensitive  than  others,  dependent  upon  the  precise  na¬ 
ture  of  the  irregularity.  A  smaller  noise  to  signal  ratio 
of  10  ~4  will  be  our  standard  value  for  the  calculations 
below;  it  describes  the  measurement  noise  relative  to 
the  power  in  the  Fourier  coefficients,  a  quantity  which 
for  white  noise  decreases  with  the  number  of  Fourier 
coefficients — by  Parseval’s  theorem.  The  value  is  ar¬ 
bitrary,  but  is  chosen  so  as  to  strongly  influence  the 
weakest  components  of  the  available  information. 

4.  Two  dimensions 

Consider  a  regular  two-dimensional  grid.  A  complete 
set  of  basis  functions  is  now  composed  of  functions  of 
the  form 


Fig.  6.  (a)  Diagonal  component  of  E  matrix,  for  one  irregular 
point,  cosine  components  only,  compared  to  that  of  completely  reg¬ 
ular  conventional  case.  The  presence  of  one  irregularly  spaced  point 
greatly  increases  the  noise  susceptibility.  ( b )  One  row  of  E  matrix 
for  conventional  aliased  case,  showing  the  alias  as  a  strong  covariance 
between  the  error  in  a  particular  frequency  component  and  the  cor¬ 
responding  alias  frequency. 


where  Lq ,  Lr  are  the  two  appropriate  length  scales,  q , 
r  are  the  corresponding  physical  positions  and  m.  n 
are  integers.  The  two-dimensional  sampling  theorem 
of  Petersen  and  Middleton  ( 1962)  shows  that  wave- 
numbers  lying  outside  the  baseband  |w,|  <  (ir/Aq), 
|  u)r\  <  (rr/A r)  will  be  aliased  into  it. 

It  is  possible  to  repeat  essentially  all  the  experiments 
done  above  in  one  dimension.  It  might  appear  that  the 
two-dimensional  case  is  of  great  interest  for  the  orbital 
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sampling  issue  as  it  relates  to  study  of  the  mean  sea 
surface.  But  consider  that  the  mean  surface  is  often 
computed  by  first  forming  a  set  of  mean  arcs,  each  as 
an  average  over  many  repeat  cycles.  In  a  1 7  day  repeat, 
such  as  that  for  GEOSAT,  any  motions  with  time  srales 
shorter  than  the  Nyquist  period  of  34  days  will  alias 
directly  into  apparent  periods  longer  than  34  days,  in¬ 
cluding  zero  frequency,  unless  removed  by  filtering  first. 
The  potential  for  distorting  estimates  of  both  the  mean 
and  the  low  frequency  variability  is  evident,  especially 
in  regions  of  strong  western  boundary  currents  with 
their  intense  short-time  scale  variability. 

We  will  thus  turn  directly  to  the  full  three-dimen¬ 
sional  case. 


a 


LONGITUDE 


5.  Three  dimensions 


a.  Real  GEOSAT 


C  LATITUDE 


Fig.  7.  (a)  Grid  pattern  in  part  of  the  western  North  Atlantic 
showing  the  distribution  of  data  obtained  during  repeat  cycle  10  of 
GEOSAT.  The  northwest  comer  is  Wank  owing  to  the  presence  of 
the  eastern  United  States  and  Canada,  and  some  of  the  descending 
arcs  had  no  useful  data.  Data  were  averaged  along- track  to  intervals 
of  105  km.  (b)  Time /longitude  sampling  for  the  data  distribution 
in  Fig.  8a.  Rapid  bursts  are  obtained  in  single  arcs,  (c)  Time/latitude 
sampling  for  Pane)  (a). 


Consider  the  data  distribution  displayed  in  Fig.  7a 
from  one  of  the  repeat  cycles  of  GEOSAT  in  the  North 
Atlantic.  The  data  are  laid  down  in  a  complex  space 
time  structure.  To  gain  some  insight  into  the  problem, 
we  show  in  Fig.  7b  the  longitude  structure  versus  time, 
and  in  Fig.  7c  the  latitude  structure  against  time.  In 
this  not  untypical  region,  the  data  arrive  in  bursts 


u> 


Ul 

z 

V) 


Fig.  8.  Contours  of  the  rows  of  the  full  covariance  matrix  for  the 
sampling  distribution  of  Fig.  7  and  the  assumptions  given  in  the  text. 
Only  values  with  magnitudes  lying  between  0.3  and  1 .0  are  contoured. 
No  negative  values  less  than  -0.2  occur.  An  element  lying  along  the 
diagonal  has  a  significant  correlation  of  the  error  of  its  estimate  with 
that  of  another  Fourier  component  if  a  contour  appears  elsewhere 
in  its  row.  A  full  display  and  discussion  of  this  matrix  is  beyond  the 
scope  of  the  present  paper;  it  does  show  that  the  major  issue  is  not 
aliasing  (at  least  within  the  restricted  wavenumber/frequencv  band 
used  here)  but  rather  uncertainty  owing  to  noise  and  the  irregular 
sampling  distribution. 


COSINE  SINE 


JOURNAL  OF  ATMOSPHERIC  AND  OCEANIC  TECHNOLOGY 


Volume  6 


lj«  -27T/4.25DAYS 


■j 


longitudinal  WAVE  number  (PER  DEGREE) 


(C) 

FKj.  9.  (a)  Estimate  of  the  variance  of  the  cosine  coefficients  of  the  three-dimensional  Fourier 
transform  of  data  distributed  as  in  Fig,  7,  at  frequency  <*>  =  -2*7(4.25  day).  Shown  as  radians/0  of 
latitude  against  radians/0  of  longitude.  Values  are  fraction  of  solution  variance;  thus  formal  error  of 
any  cosine  coefficient  would  be  ±VA’- value  shown.  These  figures  are  subregions  of  the  diagonal  of 
Fig.  8  at  fixed  frequency.  Panels  (b)-(f )  as  in  (a)  but  for  frequency  shown. 


latitudinal  wave  number  (PER  DEGREE)  LATITUDINAL  WAVE  NUMBER  (PER  DEGREE)  LATiTuDINAi  wave  NUMBER  (PER  DEGREE) 


Fig.  9.  ( Continued ) 
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closely  spaced  in  time  with  a  sequence  of  quasi-regular 
gaps  built  around  two  basic  intervals:  the  orbital  period, 
which  gives  the  time  to  return  to  the  area  in  the  next 
orbit,  and  the  longer  interval  when  complete  orbits  lie 
outside  the  region  under  study.  (The  case  shown  is  a 
real  one,  and  some  data  are  missing  and  some  of  the 
region  is  land-covered.) 

The  closely  spaced  bursts  have  a  fundamental  sam¬ 
pling  interval  of  one  second,  the  entire  along-track  burst 
occupying  0(  1  minute)  for  regions  the  size  of  that 


/—  2r  2r  2r  \ 


s2  =  cos 


—2r  2  x  2  r 

lFX  +  lry~8Jdl 


)• 


and  their  sum  (lower)  as  they  would  be  seen  in  time  if  actually  sampled 
as  shown  in  Fig.  7.  The  distribution  of  data  is  very  far  from  the 
regular  one  that  normally  leads  to  a  robust  separation  of  two  such 
distinct  cosines  through  conventional  Fourier  fits.  As  seen  in  the  text, 
even  though  separation  still  remains  a  theoretical  possibility,  the  am- 
pfitude  estimates  are  noise  sensitive. 


Fig.  1 1.  Singular  values  (squared)  underlying  the  solution  depicted 
in  Fig.  9;  there  are  a  significant  number  of  very  small  values,  whose 
effects  are  suppressed  by  the  a  priori  noise  estimate. 


shown  in  Fig.  7a  (the  spacecraft  moves  at  roughly  7 
km  s-1).  These  closely  spaced  bursts  will  behave  like 
the  closely  spaced,  irregular  data  discussed  in  section 
2,  and  one  can  anticipate  that  use  of  these  data  to  de¬ 
termine  many  of  the  time  dependent  components  will 
be  very  sensitive  to  slight  noise  in  the  observations. 

The  region  depicted  in  Fig.  7  occupies  an  area  ap¬ 
proximately  36°  of  longitude  by  15°  of  latitude.  To 
explore  the  sampling  properties,  we  sought  a  Fourier 
fit:  k  =  (-2rr/6°)  to  (2ir/6°)  in  steps  of  (2*736°),  / 
=  0  to  (2x/6°)  in  steps  of(2*7 15°),  at  frequencies  0, 
±(2*74.25),  ±(2*75.7)  ±  (2*78.5),  ±(2*717.0) 
days.  The  spatial  resolution  is  dictated  by  the  zonal 
and  meridional  dimensions  of  the  region.  This  region 
of  frequency/wavenumber  space  is  far  smaller  than 
can  hypothetically  be  fit  given  a  basic  sampling  interval 
of  7  km  along  track.  The  additional  portion  of  the  fre¬ 
quency/wavenumber  space  which  is  hypothetically 
determinable  can  be  thought  of  as  being  suppressed  by 
a  priori  assumption,  and  the  analysis  that  follows  must 
be  thought  of,  for  that  reason,  as  somewhat  optimistic. 
(If  the  suppressed  region  has  nonzero  amplitudes  the 
uncertainty  will  be  larger  than  we  are  going  to  find.) 
To  make  the  volume  of  numbers  manageable,  the  along 
track  data  was  averaged  over  a  distance  of  approxi¬ 
mately  2 10  km  to  give  a  distribution  with  half  the  data 
of  Fig.  7a,  leaving  249  data  points  in  about  650  un¬ 
known  Fourier  coefficients  (a  number  of  indeterminate 
Fourier  coefficients  were  carried  along  for  bookkeeping 
simplicity).  This  formal  underdeterminism  could  be 
removed  by  reducing  the  averaging  interval,  thus 
yielding  more  observation  points  than  the  Fourier  fit 
unknowns;  but  for  the  reasons  already  alluded  to, 
nothing  much  is  gained  by  doing  so,  the  extra  points 
being  very  closely  spaced,  and  hence  corresponding  to 
very  small  singular  values. 
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The  error  covariance  matrix  needs  to  be  displayed 
in  two  wavenumber  dimensions  and  one  frequency  di¬ 
mension,  which  is  difficult.  In  the  present  case,  there 
are  few  major  sidelobes  to  the  covariance  (little  cor¬ 
relation  of  errors  in  the  different  coefficients)  and  the 
major  part  of  the  results  can  be  obtained  from  the  di¬ 
agonal  elements  alone.  (But  there  are  some  significant 
sidelobes  and  for  complete  understanding,  one  needs 
to  examine  the  complete  702  by  702  resolution  matrix. 
An  attempt  at  displaying  it  is  found  in  Fig.  8.)  Figure 
9  displays  contours  of  the  solution  variance  at  several 
of  the  frequencies  listed  above,  for  the  cosine  coeffi¬ 
cients,  as  a  function  of  longitudinal  and  latitudinal 
wavenumber  (radians /degree).  The  convention  is  such 
that  unit  variance  means  that  the  resulting  estimate  is 
completely  indeterminate,  and  0  that  it  is  fully  deter¬ 
mined  with  no  uncertainty  (this  latter  limit  is  impos¬ 
sible  to  reach  in  the  presence  of  measurement  noise). 

If  a  variance  exceeding  0.5,  that  is  50%  of  the  solution 
variance,  is  arbitrarily  defined  as  beyond  the  limit  of 
useful  information,  certain  conclusions  can  be  drawn. 
Consider  Fig.  9e,  the  variance  of  the  Fourier  coefficients 
at  0  frequency. 

The  spatial  mean  [(0,0)  wavenumber],  and  the 
highest  wavenumbers  are  the  only  ones  which  can  be 
determined  with  any  confidence  under  the  assumptions 
we  have  made.  Because  of  the  irregular  data  distribu¬ 
tion,  and  the  assumption  of  noise  in  the  observations, 
the  system  must  rely  on  comparatively  poorly  deter¬ 
mined  data  differences  to  distinguish  these  wavenum¬ 
bers  at  zero  frequency  from  those  at  other  frequencies, 
making  the  result  statistically  uncertain.  Figure  10  dis¬ 
plays  the  way  two  pure  cosines  would  be  sampled  in 
practice.  The  extremely  irregular  way  in  which  these 
two  very  regular  waves  are  sampled  is  readily  apparent, 
and  suggests  the  difficulty  the  Fourier  expansion  has 
in  separating  them  from  each  other  and  all  other  waves 
present.  The  noise  level  can  be  reduced  somewhat  by 
using  a  second  17-day  repeat,  but  the  gain  would  not 
necessarily  be  a  factor  of  2  reduction  in  variance:  the 
wavenumber  domain  would  not  increase,  but  the 
number  of  frequencies  requiring  resolution  would 
double  from  the  initial  one  cycle  per  1 7  days,  to  one 
cycle  per  34  days,  and  thereby  require  precisely  twice 
as  many  Fourier  coefficients  to  calculate  as  well. 

The  results  just  described  are  insensitive  to  the  as¬ 
sumption  about  the  signal  to  noise  ratio;  when  the  noise 
variance  was  raised  to  10%  of  the  Fourier  coefficient 
variance,  the  diagonals  of  the  error  covunance  changed 
by  only  a  few  percent.  The  square  of  the  singular  values 
of  the  coefficient  matrix  A  are  displayed  in  Fig.  1 1; 
there  is  a  long  tail  of  small  singular  values.  Evidently 
the  0.1%  noise  variance  is  sufficient  to  suppress  most 
of  these  small  values  and  hence  the  major  part  of  the 
solution  variance.  Owing  to  the  small  slope,  the  result 
will  be  insensitive  to  the  precise  noise  variance.  Most 
of  the  solution  uncertainty  derives  from  the  “missing” 


data  owing  to  the  irregular  distribution  and  the  noise 
is  only  a  secondary  problem. 

Some  physical  insight  into  the  results  can  be  gained 
by  considering  a  plane  wave  of  wavenumber  (k,  l)  as 
it  would  be  sampled  along  a  single,  repeating  arc.  Sup¬ 
pose  the  wave  period  were  exactly  1 7  days,  so  that  al¬ 
iasing  in  time  would  make  a  sampled  wave  of  appar¬ 
ently  infinite  period.  If  ( k,  l)  is  such  that  there  is  spatial 
structure  along  the  arc,  then  the  Fourier  analysis  will 
“assign”  the  wave  to  the  band  of  wavenumbers  having 
the  same  projection  onto  the  arc.  So  for  example,  if 
the  wavenumber  vector  is  orthogonal  to  the  arc,  the 
wave  will  appear  to  have  both  zero  frequency  and  zero 
wavenumber  and  the  along-arc  spatial  averaging  cannot 
remove  the  alias  into  the  overall  (space  and  time) 
mean.  Thus  one  anticipates  a  tendency  for  waves  with 
crests  oriented  along  the  ascending  or  descending  arcs 
to  alias  into  zero  wavenumber.  The  story,  of  course,  is 
not  this  simple,  because  we  do  have  a  sequence  of  arcs, 
spaced  out  in  time  over  the  17  days  both  ascending 
and  descending,  and  a  1 7  day  wave  need  not  necessarily 
alias  at  all.  In  general,  however,  waves  with  orientation 
perpendicular  to  the  ascending  or  descending  arcs  def¬ 
initely  tend  to  be  less  well  determined  than  those  tend¬ 
ing  to  lie  along  them.  The  situation  is  potentially  very 
serious  when  the  usual  simple  averages  along  a  single 
arc  are  used  to  generate  a  “mean”  sea  surface. 

b.  Hypothetical  1 7-day  repeat 

The  results  of  the  last  sections  represent  a  real  case, 
the  actual  coverage  from  a  real  mission,  involving  ir¬ 
regular  coverage  not  only  because  of  dropped  data,  but 
because  of  the  region  of  interest  is  partly  land.  For 
studying  the  open  ocean,  one  might  prefer  to  see  the 
results  for  a  more  regular  area,  and  under  the  optimistic 
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FIG.  12.  Data  distribution  from  a  hypothetical  perfect  GEOSAT 
in  a  region  of  no  land  coverage.  Data  distribution  points  are  shown 
at  a  time  interval  of  0.055  X  orbit  period,  where  the  orbit  period  is 
about  101  minutes. 
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CU  =  -27T/4.25  DAYS 


CU=-27T/I7  DAYS 


(a) 


U=0  (THE  MEAN) 


(c) 

Fig .  1 3.  Same  as  Fig.  9,  except  for  the  hypothetical  data  distribution  of  Fig.  1 2.  Results  are  very 
similar  to  those  for  Fig.  9,  except  for  an  overall  small  reduction  in  variance  owing  to  the  presence  of 
more  data  points  with  perfect  coverage. 
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Fkj.  1 3.  (Continued) 


assumption  that  no  data  are  missing.  We  start  with 
GEOSAT  in  a  wholly  oceanic  region  with  the  coverage 
computed  hypothetically,  from  knowledge  of  the 
GEOSAT  orbit  alone,  as  depicted  in  Fig.  12.  The  same 
104  to  1  signal  to  noise  ratio  is  assumed  to  apply  and 
the  samples  were  placed  at  intervals  of  0.055  orbital 
periods  as  depicted. 

Some  of  the  result  is  shown  in  Fig.  1 3  for  a  subset 
of  all  the  periods  which  were  determined.  The  quali¬ 
tative  features  seen  for  the  real  GEOSAT  data  remain 
unchanged;  the  overall  variance  is  somewhat  reduced, 
principally  because  the  absence  of  land  in  the  region, 
and  the  assumption  of  good  data  everywhere,  means 
that  more  data  points  were  available  for  determining 
the  same  number  of  Fourier  coefficients. 

c.  Hypothetical  10-day  repeat 

The  planning  orbit  for  the  TOPEX  /  POSEIDON 
mission  is  a  10-day  repeat  with  an  orbital  inclination 
of  63. 1  °.  The  10-day  period  was  adopted  as  a  compro¬ 
mise  to  minimize  temporal  aliasing  while  at  the  same 
time  producing  a  dense  enough  ground-track  to  observe 
most  of  the  spatial  structure  of  the  large-scale  ocean 
surface  topography.  What  are  the  sampling  character¬ 
istics? 

The  10-day  repeat  occurs  in  about  130  orbits,  and 
thus  during  one  repeat  cycle,  fewer  observation  points 
are  laid  down  within  the  given  area  than  in  the  1 7-day 
repeat  which  requires  about  245  orbits.  If  the  number 
of  Fourier  coefficients  were  held  fixed,  there  would  be 
an  automatic  relative  increase  in  solution  variance  ow¬ 
ing  simply  to  the  greater  underdeterminism.  Thus  to 
obtain  a  direct  comparison  of  the  utility  of  the  10-day 
repeat  with  the  17-day  repeat,  the  10-day  repeat  was 


continued  over  the  same  number  of  orbits  (245 )  as  the 
1 7-day  repeat,  or  approximately  1 8  days.  The  data  dis¬ 
tribution  is  as  shown  in  Fig.  1 4.  The  frequency-wave¬ 
number  domain  was  kept  the  same  as  in  Fig.  1 3,  and 
the  new  result  is  displayed  in  Fig.  15. 

Generally,  the  overall  variance  is  reduced — i.e.,  the 
more  frequent  temporal  coverage  more  than  compen¬ 
sates  for  the  less  dense  surface  coverage.  The  spatial 
domain,  regarded  as  statistically  homogeneous,  is  well 
covered  in  either  case — because  of  the  diamond-pattern 
of  ascending  or  descending  passes,  all  longitudinal  and 
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Fig.  14.  Data  coverage  in  a  hypothetical  TOPEX /POSEIDON 
orbit,  where  all  data  are  obtained,  and  there  is  no  land  coverage. 
Assumed  period  was  10  days,  and  the  coverage  shown  is  that  after 
1 8  days,  to  make  the  number  of  data  points  approximately  the  same 
as  in  Figs.  12  and  13.  Again,  time  interval  is  0.055  X  orbit  period, 
where  the  orbit  period  is  about  1 12  minutes. 
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Flo.  1 5.  Cosine  coefficient  error  variance  for  hypothetical  TOPEX/ POSEIDON  l O-day  orbit  cov¬ 
erage  as  shown  in  Fig.  14.  The  more  frequent  temporal  coverage  greatly  reduces  the  unusable  regions 
(defined  as  0.5  times  the  solution  variance)  considerably  compared  to  the  1 7-day  repeat. 
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latitudinal  scales  are  observed,  and  the  denser  spatial 
coverage  from  the  17-day  repeat  is  inadequate  com¬ 
pensation  for  the  lost  time  coverage.  Thus  the  10-day 
repeat  produces  more  stable  estimates. 

6.  Final  remarks 

A  major  problem  is  the  great  difficulty  in  displaying 
the  complex  interdependencies  of  the  uncertainties  of 
frequency/wavenumber  composition  in  three  dimen¬ 
sions.  The  information  contained  in  the  full  covariance 
matrices  (or  the  separate  resolution  and  covariance 
matrices  of  the  SVD)  needs  to  be  used  in  any  quan¬ 
titative  application  of  satellite  data.  In  particular,  failure 
at  least  to  be  aware  of  aliases  in  a  specific  analysis  can 
lead  to  grossly  misleading  results  and  should  be  a  source 
of  concern  to  anyone  using  such  data. 

The  estimation  procedure  employed  here  can  be 
modified  simply  to  take,  account  of  the  actual  non  white 
nature  of  the  observational  noise — particularly  the 
strong  correlations  introduced  by  the  orbital  error. 
But  because  there  is  reasonable  expectation  that  for 
TOPEX/  POSEIDON  these  errors  will  be  reduced  to 
a  fraction  of  what  they  are  today  for  GEOSAT,  we 
postpone  this  extra  analysis  until  such  time  as  a  clearer 
specification  of  the  true  error  can  be  made. 

Within  its  restricted  terms  of  reference,  this  study 
suggests  that  the  10-day  repeat  is  preferable  to  a  17- 
day  repeat  because  of  the  reduced  solution  error  vari¬ 
ances  associated  with  the  former  orbit.  But  the  main 
purpose  of  this  paper  is  less  to  draw  a  firm  conclusion 
about  spacecraft  mission  design  than  to  set  a  context 
for  the  discussion. 

As  with  classical  sampling  theorems,  the  assumption 
was  that  all  frequencies  and  wavenumbers  are,  a  priori, 


equally  likely  to  exist.  If  one  can  demonstrate,  or  is 
willing  to  assume,  that  entire  bands  in  frequency/ 
wavenumber  space  are  empty  or  nearly  so,  the  con¬ 
clusions  drawn  here  could  be  totally  reversed. 

A  number  of  studies  are  underway  for  future  alti- 
metric  mission,  using  data  assimilation  into  dynamical 
models  (but  it  should  be  reemphasized  that  the  sam¬ 
pling  issues  raised  here  would  apply  to  any  measure¬ 
ment  made  from  an  orbiting  spacecraft).  For  any  par¬ 
ticular  model,  it  is  possible  that  a  long  repeat  cycle 
with  its  denser  spatial  coverage  does  a  better  job  of 
keeping  the  model  computation  closer  to  the  “truth.” 
One  would  have  to  presume  that  in  such  a  model  many 
of  the  frequencies  and  wavenumbers  which  cause  in¬ 
determinacy  problems  in  the  present  study  are  absent. 
They  could  be  absent  for  a  number  of  reasons — because 
in  a  nearly  linear  model  the  forcing  fields  do  not  contain 
those  bands,  or  in  a  nonlinear  model  because  the  non- 
linearities  do  not  fill  those  bands,  or  both.  As  long  as 
one  is  confident  that  the  ocean  has  the  same  spectral 
properties  as  the  model,  then  experiment  design  based 
on  the  mode!  would  be  most  appropriate.  But  the  dan¬ 
gers  inherent  in  such  a  strategy  will  be  obvious.  Space¬ 
craft  are  global  measuring  devices,  and  demonstrating 
frequency/wavenumber  consistency  of  a  model  with 
the  global  ocean  would  be  a  formidable  and  probably 
impossible  task  at  the  present  time,  given  what  is 
known  about  the  world  ocean,  and  the  limitations  on 
global  models.  With  any  new  instrument,  the  normal 
sampling  strategy  is  a  conservative  one — that  is  it  errs 
on  the  side  of  oversampling — because  sampling  theo¬ 
rems  tend  to  be  rather  unforgiving.  Determining 
whether  the  1 0-day  orbit  is  really  a  more  conservative 
choice  than  a  longer  period  repeat  will  require  more 
elaborate  studies  than  the  one  produced  here  (if  nothing 
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else,  much  greater  computing  power  must  be  brought 
to  bear  on  the  problem).  But  the  language  of  the  dis¬ 
cussion  must  be  that  of  sampling  theory. 
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